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FOREWORD 


The  long-range  propagation  of  acoustic  transients  produced  experimentally 
from  absorption  of  laser  energy  by  water  is  examined  theoretically'.  The  linear 
theory,  which  it  is  argued,  applies  for  ranges  r  >  r^  where  viscosity  effects  will 
dominate,  results  in  a  diffusion  equation  for  the  transient  pulse  form.  The 
solutions  exhibit  "solitonic"  behavior,  the  linear  1/r  peak  pressure  falloff  law 
going  over  eventually  to  a  1/r^  law  rather  than  the  much  more  severe  exponential 
law  which  applies  to  pure  tone  (harmonic)  waves.  The  scale  distance  for  the 
cross-over  to  r"2  behavior  (fresh  water),  rj  =  0.35  km»(2t,)^,  where  2t^  is  the 
transient  pulse  width  (one-half  the  period)  in  usee,  is  predicted  to  be  14  meters 
for  CC>2  laser  induced  pulses  observed  below  threshold  for  vaporization.  Above  the 
vaporization  threshold  by  a  modest  amount,  2t^  =  lpsec  and  r^  =  330  meters; 
calculated  pressure  levels  at  1.5  km  in  the  r~2  region,  are  40  db  re:  1  pbar  and 
the  pulse  width  has  spread  to  2.5  psec.  Necessary  generalizations  are  described 
briefly  for  taking  account  of  salinity  in  sea  water,  and  of  non-linear  acoustic 
effects. 
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CHAPTER  1 
INTRODUCTION 


Laser  energy  absorption  in  water  results  in  the  production  of  high  frequency- 
pressure  transients.  This  was  observed  experimentally,  long  ago,l>2  anj  { s 
expected  theoretically'.  The  mechanism  of  coupling  depends  upon  the  laser  wave¬ 
length,  and  the  light  intensity  in  the  absorption  region.  CO?-laser  light, 

\  =  10.6  um,  is  absorbed  in  a  very  thin  surface  layer  whose  thickness  scale  at  the 
lower  end,  set  by  Beer's  law,  is  of  order  10  pm.  Illuminated  spot  sizes  are  set 
by  beam  quality  and  focussing  geometry.  Assuming  direct  thermal  absorption  by  the 
water,  surface  vaporization  from  20°C  requires  a  local  energy  density 
up  =  2.65  KJ/cm-3  (=629  cal/g)  .  This  corresponds  to  a  laser  energy  threshold  of 
order  2.6  j  in  order  to  vaporize  a  10  pm  surface  layer  of  1  cm2  cross  section. 

In  recent  experiments  at  NSlv'C,  using  a  10  joule  (rated)  TEA  multimode  C0?  laser, 
varying  conditions  have  been  examined,  for  energies  ranging  from  about  a  tenth  of 
a  joule  up  to  6  joules,  and  for  spot  sizes  from  0.1  cm2  to  a  few  cm2.  These  input 
conditions  easily  straddle  the  threshold  energy  density  up,  and  the  induced 
pressure  pulse- form  characteristics  reflect  a  threshold  effect.  Higher  frequency 
transients,  100  -  400  ns  oscillations,  have  been  observed  for  lower  laser 
intensities  while  lower  frequency  pulse-forms,  2  ps  oscillations,  are  produced 
from  high  laser  intensities.  The  latter  have  somewhat  larger  peak  pressures,  and 
for  intermediate  laser  intensities  both  kinds  of  pulse  are  often  present  together; 
these  pulse-forms  frequently  show  a  high  frequency,  relatively  low  amplitude 
"leader  oscillation,"  followed  by  a  somewhat  stronger  low  frequency  main  pulse 
osc i 1 lat ion. 

Observed  pulse- forms  are  broad-band,  however,  since  only  a  single  pressure 
oscillation  cycle,  or  two,  of  the  main  pulse  bodies  occur.  (See  Figure  1).  Even 
at  rather  close  ranges  the  peak  pressure  amplitudes  are  generally  not  too  high 
for  the  application  of  linear  acoustics  (c.g.,  25  bar  at  a  distance  2  cm  vertically 
below  the  laser  spot,  in  one  instance.  See  below.  Section  2,  for  acoustic 
approximation  criteria).  At  certain  longer  ranges,  r  >  tq  say,  where  r  is  the 
distance  of  a  field  observation  point  from  the  illuminated  spot,  and  r  is  a  scale- 
length  for  assumed  dominance  of  linear  absorption  over  nonlinear  effects  (see  below), 
approximately  linear  propagation  and  absorption  of  the  pulses  is  expected.  In 
this  report,  1  present  the  application  of  linear  theory  to  pulse-form,  or,  in 

Carome,  1..  I  .,  Clark,  N .  A.  and  Moeller,  C.  E . ,  Applied  Physics  Letter,  4, 

1964,  p.  95.  ~ . 

) 

“Carome,  1..  I..  Moeller,  C.  E.  and  Clark,  A.  N.,  Applied  Physics,  40, 

19(>9,  p.  1  1(>2 . 
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other  words,  wave  packet,  propagation  and  absorption  at  long  ranges.  The 
predictions  of  linear  theory  alone,  without  need  for  invoking  possible  nonlinear 
mechanisms ,  are  qualitatively  different  from  the  case  of  a  single  frequency,  or 
pure  tone  (harmonic)  pulse.  In  the  latter  case,  the  sound  wave  aplitude  is 
attenuated  exponentially,  with  an  extinction  coefficient  B  proportional  to  the 
square  of  the  frequency,  for  high  frequencies:  8  =  koj‘,  where  w  =  2nf  is  in 
radians/sec.  But  the  radial  falloff  law  for  a  wave  packet  turns  out  to  be  a  much 
less  severe,  power  law.  The  reason  is  that  the  above  harmonic  exponential 
attenuation  law  implies  a  diffusion  equation  for  wave  packet  pulse- forms,  with 
radial  distance  in  the  role  of  "time'1;  that  results  in  pressure  pulses  at  long 
ranges  which  are  convolutions  of  short  range  forms  with  the  corresponding  Green 
function.  To  determine  the  law  of  long-range  propagation  and  absorption  for  wave 
packet  pulse-forms,  it  is  necessary  therefore  merely  to  evaluate  solutions  to  that 
equat ion. 

The  simple  diffusion  equation  that  results  from  8  =  kcj“  in  linear  acoustics 
must  be  generalized  if  relaxation  effects  are  also  present,  as  in  the  case  of  sea 
water  for  frequencies  below  100  -  200  KHz,  or  so.  Calculations  from  that  important 
generalization  are  not  given  here. 

The  key  assumption  I  have  made  is  the  dominance,  asymptotically,  of  linear 
dissipation  over  nonlinear  advection.  Pulse  propagation  including  nonlinear 
advection  effects  as  well  as  viscous  damping  can  be  investigated  to  better 
approximation  with  an  acoustic  Burgers'  equation  treatment.  A  recent  study  of 
this  equation  by  (Brighton  and  Scott-3  has  yielded  an  asymptotic  form  in  agreement 
with  that  obtained  herein.  (See  Hq.  (4.8).)  However,  experiments  readily  yield 
the  pulse-forms  for  r~i-Q,  which  can  supply  initial  data  for  linear  analysis,  and 
further  theoretical  work  is  planned  to  model  close-in  regions,  as  well  as  to 
compute  effects  of  relaxation  mentioned  earlier  for  sea  water.  Finally,  of 
great  importance  also  will  be  the  very  close-in  coupling  region  of  the  first  mm, 
or  so,  because  it  is  here  that  the  dynamical  arena  for  main  pulse  formation  is 
located. 


The  organization  of  the  report  is  as  follows.  In  Chapter  2,  I  review  the 
idealizations  of  acoustic  theory  and  identify  the  scale  length  i-q  for  dominance 
of  linear  viscosity  in  pulse  propagation,  absorption  and  distortion.  In  Chapter  5, 
I  derive  the  diffusion  equation  for  the  pulse  forms.  In  Chapter  4,  I  examine 
general  properties  of  the  solution  for  the  propagating  pulses.  In  Chapter  5,  I 
work  out  a  representative  model  example  exactly  and  display  general  features 
derived  in  Chapter  4.  In  Chapter  6,  I  generalize  to  a  complete  class  of  examples, 
illustrating  again  broad  propagation  features  identified  earlier  in  Chapter  4, 
such  as  "pure  tone"  rapid  attenuation,  and  long-range  power  law  peak  pressure 
falloff  (r'l  and  r"-)  for  wave  packet  pulse  forms.  In  Chapter  ",  I  discuss 
generalizations  and  extensions  of  the  present  work,  including  the  "sol itoni c" 
framework  for  understanding  results  of  the  earlier  linear  theory  developments. 

In  Chapter  8,  I  conclude  the  report,  and  a  mathematical  appendix  is  attached  at 
the  end.  (See  Appendix  A.) 


Brighton,  D.  (!.,  and  Scott,  J.  F.,  Phil.  Trans.  Roy  Soc.  (London)  A2H2,  15,  1  i 1  ~ G . 
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C11APTLR  2 

RUVIliW  OF  ACOUSTIC  SCALLS 


I’lie  acoustic  approximation  is  that  of  linearized  hydrodynamics,  with  the 
fluid  particle  velocity  v  small  in  magnitude  compared  to  the  sound  speed,  defined 
as 


c 


(3p/3p) 


12.1) 


4 

where  it  is  assumed  that  the  fluid  particles  undergo  adiabatic  processes.  In 
equation  (2.1),  p,  p  and  s  are  particle  pressure,  density  and  specific  entropy. 
The  Nav ier-Stokes  equation  is 


>  1  >  1  2  -> 

v  •  7  I  v  =  -  —  Vp  +  vV  v  ,  (2.2) 

where  v  =  p/p  is  the  kinematic  viscosity,  with  p  ;  4/3  n  +  c,  the  dynamic 
viscosity,  and  n,  C  the  shear  and  bulk  viscosities,  respectively;  for  water* 
v  =  0.045  cnr/sec.  liquation  (2.2)  neglects  x,t  dependence  of  p  arising  from  that 
of  local  pressure  and  temperature (T)  through  p  =  p(p,T).  Finally,  the  flow  is 
assumed  compressible,  so  that  V  •  v  0,  but  irrotat  iona  1 ,  so  that  7  x  v  =  0 
implies  existence  of  a  velocity  potential  f. 


v  =  V<j)  (2.5) 

For  small  density  deviations  p  '  =  p  -  pQ,  from  the  ambient  value  .  (),  the 
pressure  difference,  p'  =  p  -  p(),  from  the  corresponding  ambient  value  p()  can  be 
computed  from  equation  (2.1)  as 

P'  =  co  p'  *  ^2.4) 


Landau,  L.  U.  and  Lifshitz,  L.  M. ,  Fluid  Mechanics,  (Reading,  MA:Addison- 
■  eslcy  (I’ergamon) ,  1959),  Chapter  S. 


★ 


At  20° C. 
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where  c^  is  the  acoustic  (linearized)  sound  speed.  liquation  (2.4)  holds  to  good 
approximation  in  the  presence  of  irreversibility  due  to  general  viscous  (and 
thermally  induced*)  damping,  except  for  the  extreme  high  frequency  case  of  acoustic 
over-damping,  which  can  be  safely  ignored.**  Note  that  p'  does  not  have  to  be 
small  compared  to  p^  in  equation  (2.4);  only  pVpq  <<:  1  is  assumed,  whence,  from 
equation  (2.4),  pVPqCq  <<  1-  For  water,  p'  <<  22.5  Kbar.  A  second  relation, 
assuming  also  p'/Pq  <<  1,  viz. 


_  p-  *  v0  V-*  » 


C.5) 


follows  from  equations  (2.2)  and  (2.3)  if  the  nonlinear  advection  term 
(v  •  .) v  is  ignored,  and  if  p^,  p  have  vanishing  gradients.  The  neglect  of  the 
viscosity  term  in  comparison  to  the  dominant  first  and  second  terms  of 
equation  (2.5)  is  justified  again  in  the  absence  of  acoustic  over-damping.** 

The  equation  of  continuity, 


TF  *  (  •'•  )  ■ 


P  V  «  v  , 


(2.6) 


again  to  the  assumptions  of  small  p'  and  v,  becomes 


(2.7) 


1  he  sound  attenuation  caused  by  thermal  conduction  can  be  included  by  replacing 

v  in  the  expression  for  8  (see  below)  by  vepy  =  v  +  kTp  c*r  (8/Pnc  )*-,  where  for 
water  (AIP  Handbook,  1972),  the  thermal  conductivity  ^ 

k  =  5.8  x  10‘*  erg/cm  -  sec  -  K  and  the  volume  coefficient  of  thermal  expansion 

2  =  2.07  x  10"^k‘l,  both  at  20°C,  and  where  c  is  the  specific  heat  at  constant 
pressure.  The  increase  this  causes  is  only  0?02°o,  i.e.,  10“^  cm-/sec  over 

■j  =  4.5  x  10“ “  cm~/sec. 


*  *  _  J 

This  neglect  corresponds  to  the  condition  3  >>  ft,  where  is  the  reduced 

wavelength  for  a  given  Fourier  component.  Using  the  form 

1  -  5  ’  1  - 1  2  0  - 

0  c0  =  -y  cQ  k  v,  the  condition  becomes  X  >>  tt  v/Cq~P0  A  tor 

v  =  0.045  cm“/sec.  For  comparison,  at  1  NDIz,  X  =  1.5  mm,  five  orders  of 
magnitude  greater  than  90 
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Combining  equations  (.2.5),  (.2.4),  (2.5),  and  (2.7),  one  gets 


{"’-ii) 


4  =  0. 


(2.8) 


The  smallness  condition  on  v  can  be  made  precise  for  solutions  <p  obeying 
equation  (2.8).  For  outgoing  spherical  waves,  r<|>  =  f(t  -  r/c^)  satisfies 
equation  (2.8),  so  that 


co  ;>t 


_  4<j)  t  4 


(2.9) 


For  a  pulse  disturbance  having  scale-length  X,  the  last  term  is  of  order  ■*/  r  in 
comparison  to  *4/)r,  and  is  negligible  for  sufficiently  large  r.  Combining 
equations  (2.5),  (2.5),  and  (2.9)  one  gets 


Vr/C0  "  P  /p0CC 


(2.10) 


and  the  condition  on  fluid  particle  speed  becomes  v/cQ  <<  1. 

Co  include  the  simplest  effects  of  damping,  which  is  the  purpose  of  the 
present  report,  consider  the  relative  magnitudes  of  the  advcction  and  viscous 
( second  and  fourth)  terms  in  equation  (2.2).  These  are  of  order 


.  v  •  7 )  v  |  ~  v  /  X  , 


nonlinear  term. 


vv/X~  , 


viscous  terms. 


Ihe  ratio  of  the  two  is 


|  non  linear  term  I  vX  n'X  n* 

[viscous  term  ]  v  v,,pc()  v,.()w 


(2.11) 


by  equation  (2.10),  where  o  =  c.. /X  is  the  frequency  corresponding  to  the  time  scale 
for  the  pulse.  The  condition  that  the  viscous  term  dominate  is  thus 


p'$2nvi  f  =  0. 285f (water ) . 


(2. 12) 
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it'  f  is  in  MHz,  then  p'  is  in  bars;  f  =  0.5  MHz  gives  p'  •'  112  mbar  for  water;  in 
this  case,  p'/c^c^-  =  0.7  x  ID'S.  Note  here  also  that  c  =  1.5  mm/gs  gives 
=  5  nun,  so  the  neglect  of  “,/r  in  equation  (2.9)  is  justified  for  r  ■  5  mm. 


The  dominant  ideal  acoustic  terms  in  equation  (2.2)  are  of  order 


jv  ~  cp ' ! .  oco  > 


first  term, 


-  -P 


p"/‘  <r  ~(JP"/‘-,oco  ’ 


third  term. 


In  relation  to  these,  the  nonlinear  advection  term  is 


,  > 

nonlinear  term  v'/'K 


acoustic  terms  |  o>v 


(2.15 


0  ^oco" 


So,  for  a  moderately  strong  p*  =  22.5  bar,  the  nonlinear  term  is  only 
(22.5/22.5  x  10''’)  =  10'  times  the  dominant  ideal  acoustic  terms. 

When  it  can  be  assumed  the  peak  pressure  disturbances  do  satisfy  the 
condition  for  the  acoustic  approximation,  then  (2.12)  identifies  a  scale  distance 
for  validity  of  a  linear  absorption  and  propagation  model,  by  the  requirement 
that  the  peak  pressure  disturbances  at  r^  satisfy 


■pkUV  =  2™*>0f 


(2.14 


where  f  is  a  suitably  chosen  frequency  scale  value  (Courier  component)  for  the 
disturbance  at  r  =  r  .  For  smaller  values  of  r  the  propagation  may  be  acoustic 
to  good  approximation,  but  the  viscous  damping  deviations  from  "ideal"  behavior 
should  be  less  then  the  nonlinear  deviations.  For  r  >  r  ,  however,  linear 
damping  should  dominate  the  nonlinear  effect  unless  a  se 1 f- susta i tied  equilibrium 
ot  viscous  and  advection  effects  conspires  to  yield  a  sustained  favlorlike  shock 

I  will  examine  the  propagation  of  pulses  from  r  =  r.  to  longer  ranges 
a-- smiling  that  only  the  viscous  effects  need  be  considered.  The  results  will  be 
seen  to  agree  with  the  asymptotic  form  found  by  Brighton  and  Scott,-'  in  their 
analysis  of  the  acoustic  Burgers  equation  (see  Section  7  herein). 
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CHAPTER  5 

LINEAR  ABSORPTION  AND  PROPAGATION  (r  >r  ) 

o 


A  dispersion  relation  can  be  derived  for  the  linearized  equations.  Lisin; 
equation  (2.4)  to  eliminate  p ' ,  equations  (2.2)  and  (2.7)  when  linearized  are 

-> 

r  —  +  — Vn'  -  vV2v  =  0 

1  t  p  „  ; 


jjf  2  ->-  =  0 

— * +  ft  c  V  *  V 

1 1  *  oo 


Assuming  a  waveform  exp  i(k-r  -  oot)  , 

,  .  p2  ->•  .  -  lp  ^  n 

(-iu)  +  vk  )v  +  ipg  kp  =  0 

2-> 

ip0c  k-v  -  xuip  =  0 


In  the  far  field,  kr  >>  1,  k  =  kr  and  v  =  vr,  where  r  is  a  unit  radial  vector 
equations  (3.5)  and  (5.4)  become 

(- iw+vk2)v +ip0*kp'  =  0  (3.5) 


1 

ip  e“kv-iwp'  =  0 


(3.6) 


which  give  the  dispersion  relation, 


-1, 


det 


-iw+vk”  ip  k 
0 


Scok 


-  10) 


=  0, 


having  the  solution 


i  _  1  /  -> 

k  =  -  ll-iWc“)  /_ 
cn  0 

2  7 

.10  u)  _  U)  l 

-  —  +  1  ■  —  V  -  —  +  1  KGJ  . 


(3.7) 


(3.8) 


(3.9) 
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This  dispersion  relation  implies  a  diffusion  equation  for  travelling  wave  pulse 
forms. 

To  derive  this,  notice  first  that  v  =  v  ,  p',p'  ,  at  large  r,  all  have  the 
form  i"  1 r f t  -  r/c^) ;  e.g. 

p'  =  ^  f  (t  -  r/cQ),  kr»l,  (5.10) 


for  the  case  of  no  damping.  (See  equations  (2.9)  and  (2.10).)  The  case  *.  -  0 
in  equation  (3.9)  gives  k  =  w/c  ,  which  implies  the  equation 


3f 

3r 


Df 

at 


(3.11) 


since  it  holds  for  every  Fourier  component;  equation  (3.11)  is  obeyed  by 
rp'  =  A  f(t  -  r/c  )  .  For  <  j-  0,  equation  (3.9)  gives 

af  i  af  a2f 

c  +  K  2’ 

ar  o  at  at 


(3. 12) 


while  the  variable  change  (r,  t)-»-(r' ,  z ') ,  where  r'  =  r  and  t'  =  t  =  t  -  r/c^, 
so  that 


3  _  a  _  i  a  _ a_  a 

3r  ~  Or'  Cq  at'  ’  at  ~  3t'  ’ 


(3.13) 


gives 


Df  _  aff 

3r  2  * 

Ox 

where  I  dropped  the  (formal)  prime  on  r  in  the  last  line. 


(3.14) 


liquation  (3.14)  is  a  diffusion  equation  for  the  acoustic  amplitude,  with  r  in 
the  role  of  the  "time"  and  with  "spatial"  variable  x  ^  t  -  r/c^  being  the  physical 
time  variable  for  the  pulse  form  itself.  Thus,  asymptotically, 

(3.15) 


where  f  satisfies  equation  (3.14).  Typically,  r  is  already  in  the  far-fiold,  so 
near-field  "corrections"  to  equation  (3.15)  are  physically  uninteresting. 


Let  f  =  f (r  ,x  )  f 
initial  value  problem  f 
equation  (5. 1 4 )  to  be 


(t)  be  the  pulse- form  at  r  =  r  ;  the  solution  to  the 
=  f  for  f(x,r),  the  pulse  form  at  r  >r(),  is  given  from 


f  ( r ,  t  )  =  ( 4ttk p) 


-1/ 


2  ft 


dx'f  (x ') exp{ - (x  '-x) ^/4kP  \ 


( 3 . 10) 
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where  .  r  - 
that  spherical 

lo  check  this, 

t 


So  equation  (3 
+ 

f  dt 


r  .  In  evaluating  equation  (3.16)  it  is  important  to  remember 
sound  waves  have  zero  impulse;1'  this  follows  from  equation  (2.5), 


(3.  17) 


we  can  substitute  equations  (3.15)  and  (3.16)  into  equation  (3.17) 


'fQ(  .t'J ( 4-cpj 


exp 


d'-- 


74. 


Ar 


d  1  "  '' ()  ( 


(3.18) 


.17)  will  be  met  if  it  is  met  at  r  ,  i.e.,  if 

'f.U  ')  =  0  (3.  10) 
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CHAPTER  4 

GENERAL  PROPERTIES  OF  THE  SOLUTION 


For  pulse  forms  fQ ( t )  having  finite  width  tq,  equation  (3.16)  implies  another 
scale  length  r^  >  r  ,  defined  by  the  condition  that  equal  the  width  of  the 
Gaussian  propagator;  which  is  2  x  (2<p)^//",  i.e.. 


ri  -  ro  -  pi  =  t8K}_1  To 


This  can  be  written  also  as 


1  0-1 

r  = 

o  "1  8  Ko 


where 


i  ^  **■  (.  •  )  , 


(4.1) 


(4.2) 


(4.5) 


i n  which  ./2"  is  the  pu 1 se  bandwidth.  The  three  spatial  regions,  divided  by 
r()  and  r^,  are  shown  in  Figure  2. 

It  v  •  r  ,  the  Gaussian  is  narrow  compared  to  fp(t)  and  behaves  like  a 
-function  in  equation  (5.1(>)  if  f^(x)  is  smooth  enough, 

i I  ' 1  ~  exp  £-(  '  -  t)"/4<o1  —  ■  »  6 ( x  -  x ' )  ,  (4.4) 

i  .  i  n  r 


IT  r 


V 


r  <<  r 


1  ’ 


(4.5) 


i  \\  M>Muli  ii:  to  no  distortion.  From  equation  (5.15), 


r,  t  1 


r 


(4.6) 
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Micro  r  >>  r^,  tlic  (iaussian  is  broad  compared  to  t'() , 
about  t ^  =  0.  Now  equation  (5.16)  gives 


and  may  be  expanded 


f(M)  =  (  4  tt  k:  o  j 


dt”  t'Q(  t')jexp(-T“/  )•<(  )-i  '  fexpl  -  i“/ 4  .  .  .  j. 

r  T  /  r  +uo  V  (4.7) 


(  4  IT  K  P  ) 


-1/2  d 


dx 


■  1  l  f 

exp(_-x“/ 4<o) 


di  '  t  f()(r 


-)• 


neglecting  higher  terms,  and  where  I  have  made  use  also  of  the  zero- impulse 
property,  equation  (5.19).  Hvaluating  equation  (4.7), 


f(r,  r) 


VT 


■) 


-5/2 


exp (- x~ / 4c, ) 


r  >>  r. 


(4.8) 


which  agrees  with 
pulse  are  at  t  = 
amplitude  is 


the  form  found  bv  Crighton  and  Scott. The  two  peaks  of  the 

i  /  o  „  i  /  *>  1 

+  (2<o)1/-,  tor  an  overall  width  of  2(2ki  while  the  peak 


f  k  =  f(r,(2<p)1/2)  =  2k 


rx  •  (4.9) 


From  equation  (5.15),  taking  o  =  r  -  =  r,  one  finds 


pk 


2A  K~'/W  (ffT'  T'  f0(T')) 


(4.10) 


so  tlie  very  long  range  peak  pressure  falloff  law  is  r  The  energy  in  the  pulse 
involves  the  integral  of  p'-  over  all  i,  and,  owing  to  the  spreading  of  the 
packet  in  equation  (4.8),  goes  like 


If  f  (  t)  were  a  pure  harmonic,  =  C  sin  <  t  <  +a~,  Au  would 

be  zero  and  equations  (4.2)  and  (4.5)  would  give"  r^  =  o,\  In  this  case, 
r  r ^  is  impossible.  Although  r  <<  r ^ ,  accordingly,  is  well  met,  the 
(Iaussian  does  not  act  like  a  5- function  and  equation  (4.6)  doesn’t  apply 
either,  liquation  (5.16)  can  be  evaluated  exactly,  and  one  gets 


Vrighton  and  Scott  15,  19.-9. 


1  1 
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f(r,x)  =  {.  4  tt  <■  p ) 


=  (_4tt<p  ) 


I/?  f  r  2  i 

p)  “  /  dt '  Csin  exp  -(t'  -  t)  /4kp  I 

—  oo 

/+CO  -1 

dt"  exp  (iw^t"  -  (4k^)  1  i"~  I 

-  o o 

r  2  i 

m  Cexp  (iw^t  -  ww^p)  f  =  expt-tf^p)  f^( r),  c’g  - 


which  is  the  correct  pure  tone  attenuation  law. 


(4.11) 


From  a  practical  standpoint,  to  produce  large  values  of  r^  in  equation 
(4.2),  which  gives  sound  levels  like  those  in  equation  (4.6)  for  longer  ranges, 
equation  (4.5)  shows  that  what  is  needed  is  small  pulse  bandwidth  Aw.  Equation 
(4.11)  then  suggests  further  that  a  "pure  tone  carrier"  on  the  pulse  will  be 
attenuated  expotential ly  at  a  rate  appropriate  to  its  frequency. 

I  will  exhibit  these  qualitative  features  explicitly  in  exactly  soluble 
cases  in  the  next  two  sections. 
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CHAPTER  5 

REPRESENTATIVE  EXAMPLE 


A  convenient  analytical  model  of  the  pulses  observed  at  long  ranges 
experimentally*  is  (C  <  0) 


2  2 

f q  1  t J  =  Cx  exp(-x  /2 1  ),  -°°  <  t  <  +  co, 


(5.  1) 


for  which  equation  (5.16)  can  be  evaluated  exactly  (see  Figure  3).  We  have 


=  14"kp) 


2CI 


di't^exp 


1  ~>  l  ? 

a  •  s  u"  •  T)‘ 

1 


j 

1  ,  > 

.  -3/2 

Ct  •  1 

1  *  2‘p 

^  *1  } 

1  exp 

y.;  (,-*)] 


(5.2) 


The  width  of  fg(t)  is  the  distance  between  peaks  at  t  =  +_  t^  ( 2 1 ^  is  half  the 

period  of  the  single-^ycle  oscillation  shown  in  Figure  1),  so  equation  (4.1) 
gives  c  ~  (8k) “ ^ (2t  ) "  =  t^“/2<,  which  conveniently  represents  (but  somewhat 
underestimates**)  tne  radial  distance  for  transition  of  the  form  in  equation  (5.2) 
between  short  and  long  range  limits.  Evidently,  equation  (5.2)  displays  the 
limiting  behavior  predicted  in  Chapter  4  for  r  <<  rj  (equation  (4.5))  and 
r  >-  rj  (equation  (4.8)).  Over  the  intermediate  ranges  the  pulse  form  (5.1)  is 
preserved  here,  and  the  pulse  spreads,  having  width 


At  shorter  ranges,  the  pulse  looks  like  a  sawtooth.  Sec  also  Reference  7. 
**See  below,  page  16  (eg.  (5.. 8)). 

'Crighton,  l>.  0.,  and  Scott,  .J.  F.,  Phil.  Trans.  Roy.  Soc .  ,  (London)  A  2  P  2 , 
18,  1 97'.).  "  ' 
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Ihc  juilse  bandwidth,  Af  =  Aw/2tt,  is  a  function  of  p. 


\  f  =  1  /  2  TT I  . 


+  2<p/t“  ^ 


typ  lc 
focus 
r  ks  0 . 
f~  ( 

To  go 
K  =  0 

Summa 

press 

where 

focus 


for  one  unfocussed  beam  condition,  2tp  =0.2  nsec,  and  peak  pressures  are 
ally  1  bar  at  rssO.lm,  vertically  beneath  the  laser  spot.  For  certain 
sed  beam  conditions,  2tp«l  psec,  while  1  bar  peak  pressures  are  seen  around 
5m,  again  vertically  beneath  the  spot.  From  equation  (2.14)  with 
4tp)~^,  these  values  correspond  to  rQ(unfoc.) ~ 0. 14m,  and  rp( foe . ) ~ 3. 5m. 
t  r^,  take  =  2t^  in  equation  (4.1)  so  that  pp  =  tp/2«,  and  use 
.6  x  10”  1 7sec2/cm;  this  gives  r .  (unfoc . )  « p ,  =  8m  and  r,  (foe. )  «  p-i  =  200m. 


.6  x  10  sec^/cm;  this  gives  r^ (unfoc .) «p  =  8m  and  r ^ ( foe . ) « p p  =  200m 
ricing,  for  the  pulse  in  equation  (5.1),  linear  theory  predicts  1/r  peak 


re  falloff  out  to  r^,  roughly,  and  l/r2  falloff  at  longer  ranges  beyond- 
rp~8m  under  unfocussed  conditions,  and  r^  '  200m  under  certain  typical 


focussed  conditions. 

Combining  equations  (5.15),  (4.1)  and  (5.2)  (with  C  =  -  1),  the  expression  of 
the  acoustic  pressure  is 


(1  +  p/p,)' 


•  exp 


1  2/^2 

-  2  T  ftl 


(1  +  P/P, 


r  ( 1  +  P/Pp.) 


-T  1 

•  —  •  exn  -.7 

r>  t  '2 


x2/tp  (1  +  p/p. 


where  p^  is  the  peak  pressure  at  r  =  r^.  The  peak  pressure  for  r  >  r^  is 

.  =  po  ro  f5>- 

p<  r(l  +  p/pp) 

hquation  (5.7)  is  plotted  in  Figure  4  for  the  focussed  beam  example  above: 

ry  =  5.5m,  p '  =  142  mbar,  =  200m.  The  cross-over  point,  of  the  intersection 

of  r"l  and  r~2  asymptotes,  is  at  r '  =  /e" #  p  =  1.649  =  530m. 

2  .... 

Since  is  proportional  to  t~,  it  can  be  extended  considerably  it  t^  is 

increased  experimentally.  Numerically, 


t —  2 

r'  =  /e  •  p ^  =  0.55  km*(2tj)  ,  t^  in  nsec. 

Conversely,  ,-j  can  be  reduced  considerably  by  reducing  tj  experimentally. 
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This  identifies  the  practical  problems  of  (i)  pulse  forming  in  the  very  close  in 
coupling;  region,  at  the  spot,  and  (ii)  what  pulse  forms  might  be  best  for 
desired  pressures  at  long  range.  In  particular  smaller  period  pulses  can  be 
studied  more  conveniently  experimentally  due  to  the  smaller  scale  required  by  p. 
Liquation  (5.7)  is  plotted  again  in  Figure  5  for  the  higher  frequency,  unfocusseu 
beam  example  above:  Tq  =  14  cm.  Pq  =7.4  mbar,  and  p^  =  8m.  In  this  way,  the 
theoretical  predictions  of  the  asymptotic  linear  theory  described  in  the  present 
report  can  be  confronted. 
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CHAPTER  6 


MORI:  GIi.MHR.4L  PULSIi  FORMS 


The  example  equation  (5.1)  is  a  special  case  of  a  general  (and  mathematically 
complete)  class,  which  I  will  display  in  this  section,  in  part  for  possible  later 
utility.  In  preparation  for  this,  consider  again  equation  (3.16),  but  in  the 
frequency  domain  now  f(p,x)  is  a  convolution,  so  its  Fourier  transform  is  a 
product.  Let 


f(p,x) 


-t-oo 

f  dm 

J  2  TT 


f(p,m)e 


stand  to  define  the  Fourier  transform  conventions,  as  well  as  f(p,w) 
equation  (3.16)  gives  immediately 


Inverti ng 


f(r,w)  =  fQ(oj)  e 


(6.2) 


which,  by  the  way,  constitutes  a  recovery  of  the  harmonic  exponential  attenuation 
law.  Combining  equations  (6.1)  and  (6.2),  we  get  another  formula  for  f(p,x), 
equivalent  to  equation  (5.16),  and  which  can  also  be  derived  directly  from  the 
harmonic  attenuation  law  by  the  superposition  principle. 


...  ,  f  dm  y  .  ,  iu>x- 

Ut-'.T)  =  /  —  fQ  (w)  e 


(6.5) 


The  zero- impulse  condition,  equation  15.19),  requires,  finally. 


1,(01  -  0  . 


lor  the  example  in  the  last  section,  equation  (5.1),  one  has 


J  *  '  Ot,  .  exp 
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I  v identic,  f  (  )  and  arc  exactly  the  same  functional  forms.  This  is  a 

M  r>  'pc  rt  v  shared  by  all  the  ci  yen- funct  ions  of  the  quantum  mechanical  harmonic 
oscillator  owing  to  the  symmetry  of  the  Hamiltonian  under  interchange  of  position 
oid  momentum,  \iz.  of  11  -  l/2(pf  +  x~)  .  The  eigen- funct i ons  are  products  of  a 
c  mui.on  (iauss  i  an  factor,  exp  i-  -  '2tT  I  in  the  present  notation,  and  Hormitc 
•  •••!'.  n«r;iia  Is ,  Hi  / 1  j  )  ,  m  =  0,  1,  ...  It  can  be  shown  that  the  zero-impulse 
■'libitum  selects  the  odd-parity  oscillator  ei yen- funct ions ,  allowing  only 
1,  3,  3,  ...  l.quations  ( 5 .  i  )  and  16. 3)  are  the  lowest,  m  =  1,  eigen- 
i met  Me.  In  fact,  any  moderately  well-behaved  function  f~(-j,  tending  to 
.  em  at  *■  •  sufficiently  rapidly,  and  satisfying  equations  [a.  Ilf),  can  be 
c  o'anded , 


f 


where,  exp li city. 


H  i  -  / 1  .  ) 

ni  I 


c 

m 


/ 


dtf()i  ,) 


'mU/tlJ 


(6.  (i) 


(6.7) 


t 


lor  s  i  a.',  i  c  component 
n  ,  v  i  .. .  for 


"lie rm i  t  e  - pu  1  se"  forms  , 


having  c  4  0, 
in  1 


f o r  till 


IM  t .  ) 

n  1 


exp 


(.6.8) 


'  in  function  f()  possesses  exactly  1/2*' n  +  1J  cycles.  Tor  example,  f  ^  (t) 

:n  equation  la. 1 i  goes  through  a  single  full  oscillation  cycle;  for  n  =  3,  the 
.  i  ci":.  itc  polynomia  1  is  an  alternating  cubic  and  f  _(.*)  exhibits  two  oscillation 
c>  c  1  es ,  the  f.aussian  factor  pulling  the  polynomial's  down  toward  zero  at  infinity 
1  lie  cases  n  =  3  and  3  are  shown  in  figure  (>.  The  n  =  3  case  is  not  a  bad 
representative  of  the  lower  curve  in  figure  1;  but  for  large  n,  the  curves  will 
t  c-iiel  towards  symmetry  about  the  mid-point,  where  the  peaks  will  be  smallest,  and 
these  are  rather  "unlikely"  individual  shapes.  Nevertheless,  the  large  n  curves 
are  pu 1  so- carriers  of  increasing  overall  (pulse)  width,  of  order  2n^'~tj  so  they 
are  useful  for  a  brief  propagation  study  here.  (The  carrier  frequency  is  of 
order  l/2*(n  +■  l)/2n^/“tj  =  fn/tt  .  ) 


The  propagated  llermi te-pu 1 se  forms,  from  equation  (6.8),  can  be  evaluated  in 
closed  form  from  equation  (6.3).  1  do  this  in  the  Appendix;  the  result,  equation 
(A. 13)  with  the  constant  c  re- inserted,  is 


see,  for  instance,  laigen  Merzhacher,  Quantum  Mechanics,  (New  Yorktdohn  Wiley  fi 
fins  Inc.,  IPO]),  pg .  36  et  seq.,  for  properties  of  the  oscillator  e i gen- funct ions 
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where 

•  (•:  •  '  )  ■ 

l.quation  1.6.D)  for  fnl,  ,t)  reduces  to  fon(,T)  for  =  0,  and  displays  l/2«(n  +  1) 
oscillation  cycles,  accordingly  .  for  n  >  1,  the  oscillations  disappear  completely 
beyond  .  =  t^/dc  since  this  is  the  value  for  which  the  common  agruroent  of  the 
square  roots  in  equation  io.‘>;  changes  sign.  At  .  =  ti-/dr,  wc  have 


where  1  have  used  iln(.sj  =  dn£n  +  0(:.n  )  for  large  ", .  Peak  values  of  this 

pulse-form  expression  occur  at  t  =  t(2n)l/’tj.  Beyond  ,=t]~/dr,  the  argument  of 
lln  is  imaginary,  and  all  the  coefficients  of  the  polynomial  take  the  same  sign 
(after  extracting  a  factor  i.)  See  equations  (A.  17). )  So,  all  that  is  left  of 
the  llermite-pulse  is  two  endpoint  t ransient-pulse  oscillations;  the  (n  -  1)  cycle 
"llermi  tc-carrier"  center  of  the  pulse  is  gone'.  See  Figure  7. 


IB 
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CENTRALIZATION  AND  FXTF.NS10NS 


The  observation  that  acoustic  pulses  undergoing  linear  attenuation  with  an 
extinction  coefficient  proportional  to  w  should  obey  a  diffusion  equation  is  not 
new,  and  was  made  long  ago  by  Landau  and  Lifshitz.-’ 

However,  in  the  literature  of  subsequent  acoustical  theory.  Burgers' 
equation , 


n.n 


has  been  used  to  represent  combined  nonl  inear ,  and  linear  viscosity  propagation 
effects.10  The  advantage  to  this  is  that  equation  (7.1)  can  be  solved  exactly.11 
For  one-dimensional  geometry,  Blackstock  has  derived* 


>u 

•x  '  C0 


2  (  ,  B  \  ait  rU 

I  • 


(7.21 


where  B/A  is  the  nonlinearity  parameter  arising  from  the  equation  of  state  p  = 

Pi  ,s)  , 


I'  ■  l’„  *  vn‘  ■-  •  *  1 


C.3) 


and  where  u  is  the  fluid  particle  velocity.  As  previously,  r  =  t  -  c  x 
is  the  retarded  time.  Notice  equation  (7.2)  reduces  to  equation 
(a.  14)  if  t he  nonlinear  term  is  dropped. 


For  spherical  waves,  a  more  complicated  form  taking  into  account  additional 
nonlinearity  effects  associated  with  spherical  spreading,  and  neglected  in  the 


‘  bee  Reference  page  30 1 . 

1  o  .  .  , 

i  righton  and  Scott  .13,  1  <_) 7 1) . 

I  ! 

(:-  1!-.  Linear  and  Nonlinear  have'.,  f\ew  York: John  Wiley  Son-.  Inc., 
I  •*’'  t )  ,  Fhapt er  1 . 

I  he  most  readable  source,  with  primary  references,  is  Hover's  "Nonlinear 
icons  t  i e s . ” 
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analysis  I  have  given,  was  also  derived  by  Yiugo  I  '  n  y  k  h  and  co-workers  "and  by 
Rlackstoek.il  Starting  from 
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(7.4) 


from  which  again  equation  (3.14)  for  f  =  ru  follows  if  the  nonlinear  term  is 
dropped,  a  modified  Burgers'  equation  was  derived  for  the  quant itv  IV  =  u  • 
log  r/ry,  where  r^  is  a  "source  radius,"  and  where  the  independent  variable  r  is 
replaced  by  a  "stretched  coordinate,"  proportional  to  log  r/rq.  In  the  modified 
form,  '  in  equation  (7.1)  is  replaced  by  a  function  of  r,  and  the  resulting 
"Burgers'  equation"  has  not  been  solved  exactly.^ 

These  generalizations  can  be  undertaken  for  the  present  problem,  but  the 
extension  to  the  more  general  form  of  3  (w) ,  to  take  account  of  the  knees 
associated  with  onset  of  chemical  relaxation,  of  MnSO  d i sassoc i a t ion  and 
recombination  primarily,  in  sea  water  is  also  of  interest.  See  Figure  8.  The 
dispersion  relation  equation  (3.9)  is  replaced  by 


1  9 

where  is  a  relaxation  time,  and  k '  is  proportional  to  the  salinity  “  S; 


AS 


s  0.092S, 


( ~.0a) 


with  \  and  H  i  no  id  constants^,  and  S  in  ppt ;  in  F.q .  (7.5)  ,  also,  P  is  the  ambient 
■  ■!v---ure  and  1.03!  \  1  O'0  (  bar  1  "  ^  .  In  addition,  is  temperature  dependent 
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to  dominate  equation  (7.: 
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(.7.7) 
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since  now  =  (f/f_)“  <<  1.  For  T  =  25°C  and  S  =  55  again, 

k" ~  25k .  S  S 

In  the  linear  regime,  equation  (6.5)  is  readily  adapted  to  sea  water 
application,  without  finding  the  differential  equation  to  which  (”.5)  cor  responds . 
In  fact,  one  has  simply,  from  linearity, 


f(p,t) 


f0tm)eiu-'T-Btw)p 


(7.8) 


where  for  B(uj)  one  inserts  the  imaginary  part  of  k  from  equation  (7.5).  For  the 
nonlinear  viscous  case,  however,  this  procedure  is  too  simple,  and  can  not  be 
generally  applicable. lb 

Returning  to  the  simple  Burgers'  equation,  it  is  fundamentally  related  to 
the  linear  diffusion  equation  because  the  equations, 


JJJ 

)  T 


uU 


(7.9) 


•HJ 
)  x 


comprise 


a  Back  land 


trans  fo  rniat  ion 


18 


imp  lying 


)  u 

)  X 


■5 


(7.10 


(7.11 


liquation  (".11)  is  Burgers'  equation,  of  course,  and  its  solutions  are  constructed 
from  those  of  equation  (".12)  In  means  of  either  of  equations  (7.9)  and  I". 10). 

^’Crighton  and  Scott  15.  19~9. 

(Vhitham,  0.  It.,  "Linear,"  Chapter  1. 

18 

Ablowitc,  Mark  .1.  and  Segur,  Harvey,  Solutions  and  the  Inverse  Scattering 
itjuisform,  (Phi  1  ade  1  ph  i  a  :  S I  AM,  1981)  ,  Chapter''.  ’ 
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['he  soiitonic  behaviour  commonly  observed  in  problems  admitting  this  kind  of 
treatment  is  exemplified  in  a  theorem  applying  to  completely  integrable  systems:^'1 
if  a  solution  of  an  initial  value  problem  on  <  x  <  +  °°  contains  "solitons"  then 
as  r  *  -  the  "solitons"  remain  0(1).  For  acoustic  transients,  f  =  rp  is  0  (r-*) 
asymptotically  which  is  in  between  0(1),  which  holds  in  the  absence  dissipation, 
and  0(e.\p  -  hr),  which  holds  for  pure -tone  pulses  undergoing  linear  dissipation. 


U’  \blowit and  Scgur,  "Solutions,"  Chapter  3,  p.b8 
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Cl  IAPT1!R  8 
CONCLUSION 


Laser  energy  absorption  in  water  produces  experimentally  observed  ultra¬ 
sonic  transients,  which  arc  predicted  from  acoustical  theory  to  exhibit 
"solitonic"  propagation  behaviour  at  long  ranges.  Beyond  a  scale  range  r^,  at 
which  viscous  effects  become  comparable  to  nonlinear  effects  in  governing  the 
propagation,  the  peak  pressure  falloff  at  first  obeys  a  1/r  law,  but  eventually, 
for  r  >>  r,  ,  a  second  scale  range,  it  goes  over  to  a  1/r-  law.  This  result 
agrees  with  that  of  a  recent  theoretical  analysis  of  the  asymptotic  solution 
to  the  acoustic  Burgers'  equation  for  spherical  geometry.-1’  Pressure  transient 
pulse  forms  obey  a  diffusion  equation  in  the  linear  regime,  with  r  in  the  role 
of  the  "time"  and  diffusion  coefficient  r  =  for  fresh  water,  where  &  is 

the  harmonic  attenuation  constant.  When  nonlinearities  are  included,  the 
diffusion  equation  gives  way  to  a  Burgers'  equation,  which  will  apply  for 
r  <  r^.  When  MgS0_j  relaxation  in  sea  water  is  important,  the  propagation  can  be 
determined  for  long  ranges,  r  >>  r(),  from  linear  theory. 

In  the  present  paper,  I  have  examined  the  propagation  of  pulse- forms  very 
roughly  representative  of  those  seen  experimentally,  and  described  their 
properties  in  detail  for  the  linear  regime,  r  >  r(),  for  fresh  water. 


''frighten  and  Scott,  la,  Id  . 


METERS 


FIGURE  4  PEAK  PRESSURE  AS  A  FUNCTION  OF  r  FROM  EQUATION  (5.7),  FOR  REPRESENTATIVE 
OBSERVED  PULSE  PARAMETER  VALUES  (DARK  CURVE)  (r,  IS  THE  SCALE  LENGTH 
FOR  CROSSOVER  tor'2  BEHAVIOR.) 


FIGURE  5  PEAK  PRESSURE  AS  A  FUNCTION  OF  r  FROM  EQUATION  (57).  USING  A  TYPICAL, 
LOWER  AMPLITUDE,  FASTER  PULSE  (DARK  CURVE)  (THE  UPPER  CURVE  SHOWS  A 
PORTION  OF  THE  PLOT  FROM  FIGURE  4.) 
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FIGURE  7.  PULSE  FORM  AT  f ■  =P,  =  (2t,)2/8x.  FOR  THE  CASE  f„  (t)  =  -  Hn  (r/t,)  EXP  -■ 
DASHED  CURVES  INDICATE  THE  ENVELOPE  OF  "CARRIER'’  OSCILLATION 
THAT  WOULD  BE  PRESENT  HAD  THEY  SURVIVED  PROPAGATION.) 


z/2tv  (THE 
PEAKS 


(^(dB/KYD) 
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f(KHz) 


FIGURE  8  SCHEMATIC  PLOT  OF  THE  SPECTRAL  SOUND  EXTINCTION  COEFFICIENT  p  AS  A 
FUNCTION  OF  FREQUENCY  f  =  co/2n  FOR  WATER  (THE  DASHED  LINE  SHOWS  THE 
FRESH  WATER  FORM  ASSUMED  IN  THIS  REPORT) 
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by  equation  (A.l),  wheioupon. 
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i-.quat  ion  ( A .  (> )  verities  the  claim  in  chapter  <> ,  that  t,,  and  its  louner  trails! 
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To  propagate  f  (.  : I  1  have  to  evaluate  equation  ((>.3), 
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i’ roc ceding,  the  right  side  of  equation  (A. 9)  is,  with  some  rc-arrangim 
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where 


h  +  2kp 


(A. 12) 


=  (ti  -^2icp^2)  v 


(A. 13) 


Combining  the  expansion  of  F(a,f)  from  equation  (A.l)  with  equations  (A. 9)  and 
(A.  11)  - (A.  13)  , 
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(A. 14) 


whence 


exp(-^x2/t2) 


(A. 15) 


for  the  general  form,  equation  (6.8),  one  has  then  also 


f(p.t)  =  L  cn^n  ( P  *T ) 

n : odd  n  n 


(A. 16) 
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where  f  (p,x)  is  given  by  equation  (.A.  15).  The  first  few  odd  order  Hermite 
nolvnonunals  are 

R 

\\  =  25  ,  H3  =  -  125  + 

l\  =  1205  -  lbOC-5  +  3255.  (A. 
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